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Mean-field model for magnetic orders in NpTGas with T=Co, Ni or Rh 

Annamaria KiSS* and Yoshio KuRAMOTO'l' 
Department of Physics, Tohoku University, Sendai, 980-8578 

66' 

' Characteristics of magnetic transitions in NpTGas with T=Co, Ni, Rh are explained in 

(N . a unified way with use of a crystaUine electric field (CEF) model of localized 5f^ electrons. 

^ • The model takes a CEF doublet and a singlet as local states, and includes dipolar and 

' quadrupolar intersite interactions in the mean-field theory. Diverse ordering phenomena are 

If^ • derived depending on the magnitude of interaction parameters, which qualitatively reproduce 

the experimentally observed magnetic behaviors in NpTGas. The quadrupole degrees of 
freedom are essential to the diverse magnetic orders. It is argued that NpRhCas is close to 
a multicritical point where quadrupoles and dipoles with different directions are competing 
to order. 
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c , 

O . 1. Introduction 
O 

Tetragonal 115 compounds with 4/ and 5/ rare earth ions are intensively studied recently 

^ because of their intriguing behavior. For example, superconductivity at Tc = 18.5 K has been 

OO , reported in PuCoGas. The Tc is the highest among the heavy-fermion superconductors. The 

^ . second highest Tc has been observed in CeCoIus at 2.3K^ with the same crystal structure. 

' In heavy fermion systems the superconductivity and magnetism can coexist, and their pos- 

'nI" ■ 

CD ' sible interplay is a longstanding problem. Recently, NpTGas systems with various transition 

00 ! 

I metal ions T attract much attention because of their diverse magnetic behavior. Among these 

systems, NpCoGas shows an antiferromagnetic (AFM) phase transition at 47K, which can 
^ ' clearly be seen as a peak in the magnetic susceptibility.^"^ It was found by neutron diffraction 

' that the ordered moments are parallel to the tetragonal c-axis, and they have the in-plane fer- 

romagnetic structure, while the interplane stacking is antiferromagnetic. Applying magnetic 
field to this system the ordered phase is suppressed and metamagnetic transition occurs at 
low temperatures. There are two magnetic transitions in both NpNiGas and NpRhGas, but 
the magnetic structures and the nature of the ordered states are very different in these two 
systems. The magnetic structure of NpRhGas is the same as that in NpCoGas below the first 
transition, but with further decrease of temperature the moments switch over to lying within 
the a6-plane by a first order transition.®"^ The discontinuous transition is shown as a jump 
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in the magnetic susceptibility and a sharp increase of the magnetic scattering intensities. 
In NpNiGas, a canted AFM structure is found at low temperatures evolving from a ferro- 
magnetically ordered phase of moments lying parallel to the c-axis.^ In this compound both 
transitions arc continuous as indicated by finite jumps in the specific heat. At the lower phase 
transition a clear anomaly can be seen in the temperature dependence of the total magnetic 
moment. Perpendicular magnetic order with moments lying within the a&-plane is realized 
in NpFeGas below the Neel temperature 118K. A further weak anomaly was found recently 
in the thermodynamic quantities such as specific heat and magnetic susceptibility at a lower 
temperature ~ 78K, which indicates a change of the magnetic structure.® The presence of 
magnetic moments at Fe sites makes the behavior more complex in this compound. 

In this paper, we explain the mechanism of diverse magnetic orderings of NpTGas systems 
on the basis of localized picture of 5/ electrons. The observed effective magnetic moment 
in the high-temperature part of the susceptibility in NpCoGas is consistent with the Np^"*" 
(5/^) configuration, but it highly deviates from the value for Np^+ (5/'^).'^ Furthermore, the 
magnetic properties of this system above the AFM transition temperature is consistent with 
a low-lying doublet-singlet crystal field level scheme, where the doublet is the ground state. 
The next CEF level above this quasi-triplet is lying at about 1200K.^ Therefore, the low- 
temperature physics should mainly be determined by the pseudo-triplet states. Due to the 
layered structure of 115-tetragonal systems, the most important interactions are within the 
a6-planes. We introduce a two-dimensional model by taking the doublet-singlet CEF model 
within the 5/'^ configuration to explain the complex properties of NpTGas with T=Co, Ni 
and Rh. The quadrupole degrees of freedom are found to be essential in understanding the 
diverse magnetic behavior of NpTGas systems. 

This paper is organized as follows. In section 2 we introduce the model Hamiltonian within 
the pseudo-triplet subspace by clarifying the relevant multipolar interactions. Three different 
limits of the model are studied in details in section 3, and properties of compounds with 
T=Co, Rh and Ni are explained in connection with these limits. The summary and discussion 
is given in the last section. 

2. Doublet-singlet CEF model 

The Hund's rule ground state of 5/^ configuration of Np^+ ions is L = 6 and S = 2, 
which gives J = 4 as the total angular momentum. In tetragonal symmetry the nine-fold 
degenerate J = 4 multiplet splits into five singlet and two doublet states. We work in the 
following doublet-singlet local Hilbert space 



\d±) = a| ± 3) + ^/l - a'^\ =F 1) 



s) = c(| + 4) + |-4)) + Vl-2c2|0), 



(1) 
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which seems to be consistent with the magnetic properties of NpCoGas in the paramagnetic 
phase. ^ We introduce A as the energy separation between the doublet ground state and the 
singlet state. The decomposition 

{E + Ai) ^{E + Ai) = 2Aig e A2u ® Big B2g ®Eg®Eu 

shows that the pseudo-triplet space carries the Jj (^2u) and [Jx, Jy] {Eu) dipoles and 
(Big), Oxy {B2g), [Ozx, Oyz] {Eg) Quadrupoles as possible local order parameters. 

The most important interactions are within a two-dimensional layer including Np ions. For 
each a6-plane, the dipole operators at site i are written as Jk^i with k = x,y, z, and quadrupole 
operators as O^^i with /x = yz, zx, xy. We consider the nearest-neighbor interactions given by 

= ~ 4 X] ( X] ^kJk,iJk,j + XI ^iJ-Oix,iOi^j I , (2) 

(ij) \ k /. / 

where {ij) is a nearest-neighbor pair within a tetragonal a6-plane. We take the mean-field 
theory where the factor 1/4 accounts for the number of nearest-neighbors. 

Taking the CEF parameters in basis (1) as a = 0.87 and c = 0.48 seems appropriate to 
describe the high-temperature magnetic properties of NpCoGas.'^ In the absence of further 
information, we use these parameters also in the cases of T=Ni and Rh. With these CEF 
parameters the eigenvalues of operators J^, Jx, Ozx and Oxy are derived as ±2.03, ±1.98, 
±2.35 and ±2.19, respectively. Since all the eigenvalues are around ±2, we normalize the 
operators and O^ to the same value for simplicity of calculation, and for transparency of 
the model. Namely, we put their eigenvalues to be ±1,0 within the pseudo-triplet subspacc so 
that the interaction parameters should roughly be multiplied by 1/4 for the estimate of their 
magnitude. 

The interaction Hamiltonian TCab leads to very complex ordering phenomena even within 
the pseudo-triplet subspace. We study the following limiting cases of the model. As the sim- 
plest limit, we include only and Axy as the intersite interactions (Case I). Then the singlet 
excited state and the doublet are decoupled. Next, we introduce nonzero dipole interaction 
Ax = Ay keeping the quadrupolar interaction Axy (Case II). Finally, we keep both dipole inter- 
actions Ax = Ay and A^, but we assume that the quadrupolar interactions Azx = Ayz dominate 
over Axy (Case III). We argue that these limits are relevant to describe qualitatively the main 
properties of compounds with Co, Rh and Ni, respectively. Table I summarizes the limiting 
cases. We consider the two-dimensional model given by (2) in Case II and Case III, while only 
in Case I we additionally introduce Ising-type interlayer coupling in order to calculate also 
properties in external magnetic field. 
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Table I. Limiting cases of the model. Dominant components of dipole and quadrupolar interactions 
are shown in each case. 





dipole 


quadrupole 


relevance 


Case I 




^xy 


NpCoGas 


Case II 


A-2, Aj; = Ky 


^xy 


NpRhGas 


Case III 




^zx ~ ^yz 


NpNiGas 



3. Limiting cases of the model 

3.1 Case I 

As the simplest limit of the model, we first consider only A^ and K^y setting A^ = A^^ = 
= = 0. Then the singlet excited state is decoupled. The relevant operators in the 
doublet are represented by the following matrices: 

J. = I J M , o., = I : „ 1 . (3) 




^ 

Thus the doublet can be diagonalized in two different ways leading to Ising-like magnetism 
{Jz) or quadrupolar order (O^j/)- The ground state is magnetically ordered with (J^) ^ for 
^z/Ky > 1) while for Az/A^y < 1 we get quadrupole-ordered ground state with {Oxy) ^ 0. If 
we neglect the interlayer interaction, the dipole transition temperature is simply given by 
Tz = Az, while the quadrupole transition temperature Tq is given by Tq = A^y. 

We argue that the situation in NpCoGa5, which has a single phase transition at T^n = 47K, 
can be described by the doublet limit with A;^ > A-r^. The energy separation A ~ 87K, as 
estimated as from the high-temperature susceptibility, should be important to the suscepti- 
bility anisotropy for example. However, the singlet does not influence the main properties of 
the model such as the pattern of ordered phases or the nature of the phase transitions. 

The AFM ordering in NpCoGas is characterized by the ordering vector K = (0,0,1/2) 
in unit of (27r/a, 27r/a, 27r/c). In the following, ordering vectors k or K always mean three- 
dimensional vectors in unit of (27r/a, 27r/a, 27r/c), while ordering vectors q or Q are two- 
dimensional vectors in unit of (27r/a, 27r/a). To discuss properties in the presence of external 
magnetic field like magnetic susceptibility or temperature-magnetic field phase diagram, we 
introduce an interlayer interaction only in the present case as follows: 

^c = ^Afj:j<;)jiy^), (4) 

I 

where we have explicitly introduced the layer index I in this case, and the factor 1/2 accounts 
for the number of nearest neighbors along the c axis. Although the Ising-type interaction has 
been chosen for simplicity, there is no difference of the result even if we choose an isotropic 
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T/A T/A 

Fig. 1. Left: Temperature dependence of the magnetic susceptibility in the Ising doublet limit. Inset 
shows the magnetization curves as a function of temperature for fields h = 0.002, 0.008, 0.017, 
0.027 and 0.034. Right: T-H phase diagram for field direction (001). Filled squares mean second 
order transitions, while empty squares represent first order one. We took Xz + = 0.586 with 
X-^/Xz = 0.05 in the calculation, which gives the zero field phase transition as T^/A = 0.54. 



interlayer interaction. The Ising-type AFM transition temperature is given by 

1-(A, + A^)xl(T,) =0, (5) 

where we have assumed > 0, and the local susceptibility Xl{T) is given by Xl{T) = 
+ exp(— /3A)/2]^^ with /3 = 1/T. The (homogeneous) magnetic susceptibility for il||(001) 
in the paramagnetic phase {T > T^) is given by 

Note that the antiferromagnetic interplane interaction A^ > reduces the homogeneous 
susceptibility. 

The left part of Fig. 1 shows the magnetic susceptibility for both magnetic field direc- 
tions (001) and (100). We use dimcnsionless coupling constants and magnetic field such as 
Afc,^ = Afc,^/A, and h = gji^H / {k^/S), where g = 3/5. If we decrease the ratio A^/A^, i. 
e., towards weaker interlayer coupling, the peak in the magnetic susceptibility for if||(001) 
becomes sharper because of the decrease of the denominator in expression (6). 

We also calculated the temperature dependence of the magnetization, which shows mcta- 
magnetic transition at high fields and low temperatures as shown in the inset of Fig. 1. The 
calculated temperature-magnetic field phase diagram can be seen in the right part of Fig. 1. 
Let us consider the mctamagnctic transition in the ground state. The staggered component of 
the order parameter {Jz)a — {Jz)b vanishes when the magnetic field exceeds a critical value 
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Her- For H < Her the ground state energy in the mean-field theory can be written as 



Within the AFM ordered phase the conditions dJ^<^/d{Jz) a = and dJ^^/d{Jz)B = gives 
{Jz)a = —{Jz)b = 1, which leads to jr< = —l/2{Xz + A^). Similarly, in the high-field phase 
we get (Jz) = 1 from (8), which gives jr> = -1/2(A-, - Aj-) - H. The AFM order vanishes 
discontinuously by increasing the magnetic field at if = He - which gives = From 
expressions (7) and (8) with ( J^) = 1 we get Her = A^. Namely H^r is given by the value of the 
interlayer coupling constant. This relation clarifies the origin of the metamagnetic transition. 

Expression (5) shows that the transition temperature Tz is related to the value of A;j + A^, 
while the critical magnetic field Her to A^. Therefore, Her can be small relative to the transition 
temperature Tz if we choose small interlayer coupling A^. Furthermore, with small A^, the 
susceptibility is sharply peaked for i^lKOOl). These features are in good agreement with the 
measured results for NpCoGas-^"^ 
3.2 Case H 

In this limit we take nonzero dipole interactions and Aj, = A^^, and assume that the 
quardupolar interaction A^y dominates over Ayz = A^x- This A^y quadrupolar interaction 
stabilizes the magnetic moment along (110) because of the third order term {Jx){Jy){Oxy) in 
the Landau free-energy expansion. Namely, the system gains the maximum energy when both 
{Jx) and {Jy) are non-zero and equal. The third order term also shows that pure magnetic 
order of the perpendicular dipoles {Jx + Jy) does not exist, since the quadrupoles {Oxy) are 
induced by the symmetry. 

The left part of Fig. 2 shows the ground state of the model with fixed \x = Xy = 2.13, 
which value is chosen as a trial, and Xzx = Xyz = 0. In the limit of Xxy, Xz » A^, the phase 
boundary between the magnetic {{Jz) 7^ 0) and the quadrupolar {{Oxy) 7^ 0) phases tends to 
Az = Xxy, recovering the doublet limit with A^ = 0. The nonzero Xx causes the deviation of the 
phase boundary from Xz = Xxy in such a way that the mixed phase with {Jx + Jy){Oxy) 7^ 
expands for Xxy < A^ at the expense of phase with {Jz) 7^ in the ground state. The right 
part of Fig. 2 shows the ground state magnetic moment 



the ordered magnetic moment decreases with increasing quadrupolar interaction, and finally 
disappears at Xxy = 3.26 in the case of A^: = 2.13. The vanishing of the moment gives the 
second order phase boundary between the pure quadrupolar phase and the mixed phase. The 



^< = \Xz {{Jz)l + {Jz)l) - lxi{J.)A{Jz)B - ^(A. + A^) i{Jz)A - {Jz)b) 



(7) 



while for H > Hnr it becomes 



^> = ^(A. - ^z){Jz)^ - (A. - Af )(J.> - H. 



(8) 




6/17 



J. Phys. Soc. Jpn. 



Full Paper 




Fig. 2. Left: The ground-state phase diagram with A^; = 2.13. For details see the text. Right: Magnetic 
moment m = {J^ + Jy)/V^ in the ground state as a function of Xxy for Aa; = 2.13 and A^ < 1.25. 

boundary is given by Xxy = 2Xx — 1, which can be seen as the horizontal straight line on the 
phase diagram shown in the left part of Fig. 2. 

The compound NpRhGa5 shows two successive magnetic transitions. From neutron diffrac- 
tion it is concluded that the magnetic structure below the first transtion = 36K is the same 
as that in NpCoGas, while at = 33K the moment changes the direction to {Jx + Jy) 7^ by 
a first order transition. The moments are ferromagnctically ordered within each 06-planc.^"^ 
The ordered moments {Jz) completely vanish below 33K. In order to simulate the situation 
in NpRhGas qualitatively, we now fix a part of parameters as Xxy = 1-5 and A^; = 2.13, and 
derive the phase diagram in the A^-T plane. In the parameter range Xxy ^ A^: the quadrupolar 
interaction prefers in-plane magnetic moments to the c-axis moments. A consequence is that 
there appears a regime on the A^-T phase diagram where an ordered phase with {Jz) 7^ 
undergoes a first order transition to another one with (J^ -|- Jy){Oxy) 7^ with decreasing 
temperature. The left part of Fig. 3 shows the calculated phase diagram. Along the path 
(A2 = 1.68) shown by the arrow, the moment /v, = ^ [Jx + Jy)"^ /2 -\- J| develops as shown 
in the right panel. The inter-plane AFM coupling A^ can simply be included in the mean- 
field theory by modifying A^ ^ A^ -|- A^. The magnetic moment increases discontinuously at 
the first order transition. This behavior is in qualitative agreement with the results found in 
NpRhGas.^ Note that with A^; -|- ~ 1.62, the three phases meet at the multicritical point. 

We note that the consideration of the singlet state is essential to obtain the perpendicular 
magnetic order in Case II. This is in contrast to the case of NpCoGas, which can be described 
in the doublet limit. However, the value of the energy separation A is not crucial to the 
existence of the phases obtained and their characteristic features. 

Although we considered a two-dimensional model taking into account only the interactions 
within the a6-planes, we comment shortly on the three-dimensional magnetic and quadrupo- 
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Fig. 3. Left: Xz-T phase diagram for X^y = 1.5 and X^ = 2.13. Dashed hnes indicate the first order 
transitions, while the sohd hne represents the second order one. Their meeting point represents a 
multicritical point which is near the arrow relevant to NpRhGas. Right: Temperature dependence 
of the magnetic moment along the arrow in the upper panel. With + = 1.68, the phase 
transitions are derived as T^/A = 1.31 and T^/A = 1.36. Inset shows the moment near the 
transitions. 

lar structure. The magnetic moments {J^ + Jy) are ferromagnetically ordered within the ab 
planes at low temperatures. This means that homogenous quadrupolar moments {Oxy) are 
present within the planes even when the quadrupolar coupling constant X^y is zero. The in- 
terplane stacking of the magnetic moments is antiferromagnetic, but the ferroquadrupolar 
order is realized on subsequent planes. This situation is consistent with the invariant term 
{Jx){'is){Jy){—iii){Oxy){0) in the Landau free energy expansion. The ferroquadrupolar order 
means that an orthorombic lattice distortion should develop below = 33K. We expect 
experimental observation of the distortion. 
3.3 Casein 

Now we assume that the quadrupolar interactions A^x = dominate over Axy in contrast 
with cases of T=Co or Rh. In this case the third order term {Jz){{Jx){Ozx) + {Jy) {Oyz)) in the 
Landau expansion plays an important role. Namely, it can lead to two successive transitions 
both of which are continuous. After the ordering of one of the multipole moments Jx, Jz or 
Ozx, a lower phase transition is possible by the simultaneous appearance of the other two 
multipole moments. A difficulty of the ordinary Landau-type expansion is that the order 
parameter of the first phase transition is not necessarily small at the second transition, so 
that the expansion of the free energy with respect to it is not justified. Therefore we use the 
method developed in ref.ll, and calculate the Ginzburg-Landau free energy functional in all 
orders of the finite order parameter which becomes nonzero below the first phase transition. 
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In the present limit wc consider the following interactions within the tetragonal a5-plane: 



ZX,lOzx,j] 



{ij) 



We start with the path-integral representation of the partition function 



11 



/ 



PX^exp 



-Sb 



/ dTHir) 
Jo 



(10) 



(11) 



where the integration variable X"^ replaces the operators J^, Jz and Ozx- Sb is the Berry phase 
term, which enters because of the non-bosonic commutation property of Jx, Jz and Ozx- We 
use the Hubbard-Stratonovich transformation for each imaginary time interval, which is the 
replacement —A^X"^ (f)^ by introducing effective fields (j)^. These fields (p"^ mediate the 
original interaction. The Hubbard-Stratonovich transformation converts the original inter-site 
interaction into a local interaction between the effective fields and multipolar operators. Using 
the property that ^'s have Gaussian distribution, the partition function can be expressed as 

r0 



-I 



where 



(12) 



(13) 



Wc further use the static approximation within which the path integral over in Z IS 
replaced by trace calculation of corresponding operators. Then we obtain the following form 
for the partition function 



Z = J P(/)^exp[-/3J^] , 



(14) 



where T is the Ginzburg-Landau free energy functional given by 



(15) 



i-,3 



with 



dl 



2 (e-^^ - 11 



d2 = -P-, 



A (2 + 0-0^)' ^ '^(2 + e-/3^)' 
and Zq = exp(— /3A) -f 2. A'^ is the number of the sites. The mode coupling free energy has 
the form 

/mc = -j^lj2 {T^ieMP^txt/2)/Zo) + \f3d,^f^f + ^/Jdi^f^f + ^/?d2</>|d • (16) 

Let us consider the case where the first phase transition at a temperature T = corre- 
sponds to the ordering of Jz dipole moments. This means that for T < Th the expectation 
value (0q) = Az{Q){Jz){Q) becomes nonzero, where A^(Q) = (l/2)A2(cos((5a:) -|- cos(Qy)) 
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is the Fourier transform of the couphng constant A^. The lower transition temperature Tl is 
derived from the condition det(x~^) = 0, where x is the generahzed susceptibihty matrix. We 
have the relation: 



X 



J-^GJ-^ - J- 



(17) 



where the interaction matrix j is given by 

/ A.(Q) 

J = 



\ 








Ax(q) 










(18) 



A^a;(-(q + Q)) / 



The matrix G satisfies the following equation 



13 



G = l3 



V 
















^q+Q'?^-(q+Q) 



(19) 



) / 



-(q+Q)'^q) (-^q 

which means the generalized susceptibility for the effective fields. We are interested in the case 
where the first transition is ferro-type with Q = 0. Therefore, the wave vector q is equivalent 
to q + Q. 

If the lower transition is also of second order, we make perturbation expansion, under 
the finite value of (^q), for the mode coupling free energy (16) in terms of order parameters 

and (f)^ up to second order. Then, the Ginzburg-Landau free energy functional can be 
written in the form 



1 1 ^ 



(20) 



where the matrix H = {Hij } composed by the coefficients of the second order terms is related 
to G as = G, and it has the same block diagonal form. We note when 0* shows Gaussian 
distribution, G^^ is given by its variance. Using the relation (17) we obtain 

A^(0)A^(q)A,^(q)i?n(i?22^^33 - H23H32) 



(21) 



[A^(0)-1 - /fii][(A^(q)-l - iJ22)(A,a.(q)-l - i?33) - i?23i^32] ' 

Therefore, the condition det(x~^) = is equivalent with det{H) = 0. Calculating the coeffi- 
cients Hij in the free energy expansion (20), we can derive the lower phase transition by the 
condition det{H) = 0. 

• T>Th 

In this temperature regime with (0q) = 0, the mode coupling free energy (16) gives 
H23 = H32 = in (20). The second order coefficient of term <^q can be read from (15) as 



1 



A, ^(2 + e-/^A) A, 



- XL : 



(22) 
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where we have used that A2(0) = A^. Due to the block form of matrices G and i7, the 
component of the generahzed susceptibihty, which corresponds to the dipole operator 
can be obtained as 

Xii = A-^(ffri^-l)= W , (23) 

which is the same as expression (6) derived in Case I with A^ = 0. This is due to the 
fact that the homogenous magnetic field lying in direction (001) is the conjugated field 
to the dipole moment Jz in the case of K = 0. The first transition temperature Th can 
be derived from Hn = 0, as the first instability in ordinary Landau mean field theory 
(see equation (21)). This leads to the condition 1 — AzXl{Th) = 0, which is the same as 
expression (5) taking A^ = 0. Or equivalently, Xu = leads to the same condition for 
the transition temperature. 

. Th > r > Tl 

In this temperature regime fmc becomes nonzero in the Ginzburg-Landau free energy 
functional (15). Under the finite value of {4>q) = Az{Jz), the Ginzburg-Landau free energy 
functional is derived as 

q 

+{A-' + d2Wo<Po-d2{rof] , (24) 



where 



(</.a) + A (</.g)-A ■ (</.g)2-A2y' 

and Zq = exp(— /3A) + 2cosh(/3(^Q)). When {^q) — we obtain 7r2 — and tti = —di, 
which means the vanishing of the mixing term as we expect. We then come back to the 
high temperature form of the GL free energy functional given by (15) with /mc = 0. 
The lower transition temperature Tl can be derived from det(x~^) = or det(^) = 0. 
Explicitly we have the condition 

[Ax(q)-' - 7ri(TL)][A,,(q)-i - 7ri(TL)] - = 0. (26) 

Similar calculation can be performed for the free energy expansion under finite {(I)q) = 
Aj;(0)(Jj;)(0) or {(f)^) = Azx(0){Ozx){0). The left panel of Fig. 4 shows the ground state 
(r = 0) phase diagram with = A^j/A fixed at 2.13. We find three regimes where pure 
dipolar (J^ or Jx) or quadrupolar (Ozx) order is realized. Thus, we can enter to the regime 
with mixed order parameters {Jz)-{Jx)-{Ozx) from these three different sides. Let us consider 
the limiting cases of the phase diagram. 
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Fig. 4. Left: Finite order parameters at T = in the Xz-Xzx plane with Aa: = 2.13. Right: Phase 
diagram in the A^-T plane with a tetracritical point. We have chosen A^ = 2.13 and Xzx = 1-5. 



(i) ^zx = 0: there is a first order transition at Xz ~ 1.25 between the two kinds of magnetic 
phases. The situation is similar to that in Fig. 2 with Xxy = 0. The phase boundary between 
the mixed and ( J^) / phases can be obtained from condition (26) by taking the limit T — 0, 
which gives 

2 







1 1 \ / 1 1 \ / 1 



1 / -^2: + 1 ~ ^zx ~ ^x 



A2 V \xXzxi\z + 1) y ' ^^^^ 

where we have used that (^q) = A;^ = A^A in the ground state. Thus, from (27) the phase 
boundary is obtained as Xzx = ^ — + ^z- The boundary between the mixed and {Jx) 7^ 
phases is derived as Xzx = — Xz{l + 2Xx)/{2Xx — 1). The intersection of these two phase 
boundaries gives the point A^; 1.25. 

(ii) A^; = 0: the point Xzx = 2.13 (= Xx) separates the (Jx) ^ dipolar and (Ozx) 7^ 
quadrupolar phases. The phase boundary between the mixed and (Ozx) 7^ phases terminates 
at 

1. 1,.. 1 



Xzx = + ^(2A, + 1) + -yJ(2Xx - 1)2 + 4A2 + 8A,A, + 12A,. 

The right panel of Fig. 4 shows the phase diagram in the A^-T plane with fixed parameters 
Xx = 2.13 and Xzx = 1-5. With A^ > 1.5, the first transition makes (J^) 7^ 0, and the second 
transition gives {Jx){Ozx) 7^ continuously. However for A^ < 1.5, the first transition gives 
{Jx) 7^ 0, and at the lower transition the term {Jz){Ozx) becomes nonzero. The transition 
temperature of the pure Jx ordering (A^ < 1.5) is given by 

A 

" ~ln(2A^ -2) -ln(2A,; + l)' ^^^^ 
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which is derived from the vanishing of the coefficient of the second order term 0f 0^ in the 
free energy expansion (15). The expression (28) leads to the horizontal line T^^/A ~ 1.8 as a 
function of with \x = 2.13 (see the phase diagram). The expression (28) also shows that 
continuous transition to the phase [Jx) 7^ is not possible at any temperature with < 1, 
i. e., < A. 

The magnetic coupling enhances the lower transition temperature by the same mech- 
anism as discussed in ref.l4. The phase diagram contains a tetracritical point at ~ 1.5, 
where the four critical lines meet. Phase diagrams with multicritical points should also have 
experimental interest, because the application of uniaxial pressure or doping can drive the 
system through a multicritical point to a different regime of the phase diagram with very 
different physical properties. 

We note that the phase diagram presented in the left part of Fig. 4 remains the same for 
the (llO)-type perpendicular dipoles instead of the (lOO)-type. Namely, there is a continuous 
degeneracy with respect to the phase transition temperature within the Jx^Jy order parameter 
space due to the tetragonal symmetry. The difference is that the ordered quadrupoles will have 
nonzero moments {Ozx + Oyz) instead of (Ozx) 7^ 0. 

In NpNiGas there appear two successive continuous phase transitions at Tc = 30K and 
Tn = 18K.^'^ The first one is the ferromagnetic ordering of Jz dipoles with the ordering vector 
K = (0,0,0). At temperature Tn non-zero (Jx) component appears with k = (1/2,1/2,1/2), 
leading to a canted antiferromagnetic structure. In ref.9 two possibilities are mentioned for 
the orderings: (i) both transitions correspond to the Np sublattice or (ii) one of the transitions 
related to the Ni sublattice. We assume that both transitions are caused by the Np ions, which 
is consistent with recent neutron diffraction results.'^ 

In order to investigate the second transition to an ordered phase with (Jz) 7^ 0, we take 
the simplifying limit Xx = Xzx- Then the free energy can be calculated by the diagonaliza- 
tion of the Hamiltonian within the basis (1). The left panel of Fig. 5 shows the specific heat 
against temperature with two second-order transitions. We also calculated the direction of 
the spontaneous moment {&) below the lower transition temperature and the total magnetic 
moment (//) as a function of temperature (see right part of Fig.5). We define the spontaneous 
moment 7^ = {Jx)i- + {Jz)^i and its deviation from the direction z as tan(^/2) = {Jx)/{Jz)- 
The magnitude of the moment is given by = {Jx)"^ + (</z)^- The calculated temperature 
dependences of 9 and jj, are shown in the right panel of Fig. 5. The lower transition causes an 
anomaly in jj, in good correspondence with the experimental result in NpNiGa5.-'^'^ 

In the case of antiferromagnetic coupling A^, < 0, the lower quadrupolax transition will lead 
to a canted AFM structure within the ab-planes. In order to have the non-zero invariant term 
{Jx){Q,x){Ozx){Q,zx){Jz){Q,z) with the ordering wave vectors in the parenthesis, the condition 
Q,zx = —{Q,x + Q,z) = ~Q,x should be satisfied with = 0. Hence the quadrupolar order 
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T/A T/A 

Fig. 5. Left: Temperature dependence of the specific heat with two transitions at T^/A « 15 and 
Tc/A Ri 13.6. The parameters are chosen as A^^ = A^^; = 16.5 and A^ = 22. Right: Temperature 
dependence of 6 (degree) and the moment /x. 

should have the same ordering vector Q^; = (1/2, 1/2) as that of Jx- 

4. Summary and Discussion 

In this work we have studied a two-dimensional mean field model composed of a non- 
Kramers doublet ground state and a singlet excited state. The model includes dipolar and 
quadrupolar degrees of freedom. The layered structure of / ions should be important in 
realizing characteristic magnetism and superconductivity in 115 systems. For theoretical con- 
sideration of these systems, a relevant picture should be a strong two-dimensional interaction 
and a weak interlayer interaction. Therefore, the main interactions which determine the mag- 
netic structure are within the two-dimensional a6-planes. The main purpose of this work is to 
understand the origin of the complex magnetic structures of NpTGas systems with T=Co, Ni 
and Rh. Thus, we have studied a two-dimensional model which includes interactions within 
the tetragonal a6-planes. We have analyzed some limiting cases of the model by changing 
the interaction parameters, and found behaviors reminiscent of the magnetic properties of 
NpTGas systems with T=Co, Ni and Rh. 

The intersite interactions in NpTGas with T=Ni may be different from those with T=Co 
or Rh since the Fermi surface of NpNiGas contains one more conduction electrons per unit cell. 
The difference in the number of cZ-electrons may be a reason why in NpNiGas the quadrupolar 
interactions Xzx = \z dominate over X^y in contrast with cases of T=Co, Rh. On the other 
hand, the CEF structure should be less sensitive to the change of the Fermi surface. It is likely 
that other Np 115 systems may correspond to a shifted parameter space. Then multicritical 
behaviors can be expected under appropriate conditions. For example, the ordering of the 
perpendicular magnetic moments can also happen first depending on the interaction param- 
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eters. The perpendicular moments experimentally found in NpPtGas^^ seem to be explained 
in this way. 

It is always a question in the case of actinidc compounds whether the electrons arc localized 
or itinerant. There is a tendency of increasing localization with increasing atomic number 
along the / series of the elements in the periodic table. The element Np, situated between 
U and Pu, seems to be at the localized-itinerant border. Recent NMR results on NpCoGas 
suggest localized behavior in contrast to the more itinerant-like UCoGas.-'^^ Among the Np 
compounds, Np02 is an example where the localized description works well to describe the 
triple-q octupolar order below 25K. Even in the case of U-based compounds, the localized 
picture may apply to the case of the famous hidden ordered phase of URu2Si2. Although our 
simple CEF model can explain the diverse behavior realized in NpTGas, the purely localized 
model cannot explain the enhanced T-linear term in the specific heat, for example. Hence 
a more sophisticated model is desirable to account for the dual character of 5/ electrons. 
We note that in the case of even number of 5/ electrons as in Np^"*", electronic states in the 
localized limit and the band limit may be connected continuously.^^ 

A further question is the value of the ordered magnetic moments at low temperatures. It 
is found experimentally that in the case of NpCoGas the saturated magnetic moment is about 
0.84/7,b/Np, and it is 0.96/iB/Np in NpRhGas. In the mean field theory, the zero temperature 
value of the moments is mainly determined by the crystal field parameters. For example, 
with crystal field parameter a = 0.87 we get 1.22/1^ for the saturated magnetic moment 
in the case of NpCoGas, which is larger than the observed value. Wc saw previously in the 
case of NpRhGas that the presence of quadriipolar interaction X^y can reduce the value of 
the magnetic moment {J^ + Jy). On the other hand, quantum fluctuations can lead to the 
reduction of the ordered moment even for localized electrons. We expect stronger fluctuations 
in the Np 115 systems because of the strong two-dimensional feature compared to the case of 
the cubic NpGas, for example, in which the ordered moments are about 1.51;Ub-^^ 

The main properties of NpTGas systems with T=Co, Ni and Rh can be understood 
within the two-dimensional model. However, in order to discuss behavior in the presence 
of magnetic field like susceptibility anisotropy, for example, the interlayer coupling should 
be also included. Only in the case of NpCoGa5 we have included an interlayer interaction, 
and calculated properties in external magnetic field. It would be interesting to incorporate 
interlayer coupling also in the cases of T=Ni and Rh and compare the obtained results to the 
measured ones. The difference in the number of conduction electrons may be a reason why the 
interlayer RKKY coupling becomes ferromagnetic in NpNiGa5, as compared with observed 
antiferromagnetic interlayer ordering in NpCoGas and NpRhGas We plan to discuss these 
issues in a subsequent work. 
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